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Abstract:  Optical  arbitrary  waveform  generation  using  the  line-by-line 
pulse  shaping  technique  has  been  shown  to  be  sensitive  to  variations  in  the 
offset  frequency  of  the  input  frequency  comb  due  to  time-domain  waveform 
interference.  Here  we  present  a  frequency-domain  model  that  is  able  to 
predict  waveform  changes  arising  from  offset  frequency  variations.  In 
experiments  we  controllably  shift  the  frequency  of  a  comb  derived  from  a 
phase-modulated  CW  laser,  which  allows  us  to  quantitatively  investigate 
waveforms  generated  by  pulse  shaping  as  a  function  of  offset  frequency. 
Experimental  data  are  in  excellent  agreement  with  the  predictions  of  our 
frequency-domain  model.  In  addition,  we  propose  and  analyze  new 
waveforms  designed  for  monitoring  of  offset  frequency  variations  by  pulse 
shaping. 
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1.  Introduction  and  problem  description 

Optical  combs  are  comprised  of  spectral  lines  that  are  evenly  spaced  and  determined  by  the 
source  repetition  frequency  and  with  well  defined  phase  relations  between  the  lines  [1], 
Optical  pulse  shaping  techniques,  in  which  phase  and  amplitude  manipulation  of  optical 
spectral  components  allow  synthesis  of  user-specified  ultrashort  pulse  fields  according  to  a 
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Fourier-Transform  relationship,  have  been  developed  and  widely  adopted  [2],  Previous  pulse 
shapers,  however,  have  generally  manipulated  groups  of  spectral  lines  rather  than  individual 
lines  due  to  difficulty  in  building  a  pulse  shaper  capable  of  resolving  each  spectral  line  at 
typical  repetition  frequencies  (<1  GHz)  for  practical  optical  combs  generated  by  mode-locked 
lasers.  This  results  in  waveform  bursts  that  are  separated  in  time  with  low  duty  factor  and  that 
are  insensitive  to  the  variation  of  the  frequency  offset  of  the  comb  lines.  However,  the  phase 
and  amplitude  of  each  individual  comb  line  should  be  independently  controlled  in  order  to 
approach  true  optical  arbitrary  waveform  generation  (O-AWG)  [3].  A  pulse  shaper  capable  of 
implementing  line-by-line  control  can  generate  shaped  waveforms  spanning  the  full  time 
period  (T)  between  the  initial  mode-locked  pulses  (100%  duty  factor).  From  a  time-domain 
perspective,  in  the  line-by-line  regime  waveform  contributions  arising  from  adjacent  pulses  in 
the  mode-locked  train  overlap  and  interfere.  The  effect  of  the  interference  depends  on  the 
relative  phase  between  adjacent  mode-locked  pulses,  which  is  proportional  to  the  comb  offset 
frequency.  Consequently,  frequency  fluctuations  of  the  mode-locked  comb  translate  into  time- 
dependent  waveform  noise  in  line-by-line  pulse  shaping.  This  effect  fundamentally  links  the 
field  of  pulse  shaping  with  the  field  of  frequency-stabilized  mode-locked  lasers  (combs). 

The  first  observation  of  such  waveform  noise  in  line-by-line  pulse  shaping  was  reported  in 
[4]  by  selecting  two  comb  lines  from  a  harmonically  mode-locked  fiber  laser  without  active 
frequency  stabilization.  The  laser  was  operated  at  ~!0  GHz  repetition  rate,  which  made  it 
possible  to  resolve  individual  lines  using  a  carefully  designed  pulse  shaper.  Selecting  only 
two  spectral  lines  from  the  broader  input  spectrum  yields  a  cosine  waveform  in  the  time- 
domain.  The  relative  delay  of  the  peaks  of  the  cosine  can  be  tuned  by  using  the  pulse  shaper  to 
control  the  phase  (O)  of  one  of  the  selected  lines.  The  pulses  are  shifted  by  T/2  (half  of  the 
pulse  period)  if  cp=7t  is  applied  to  a  single  line.  The  experiments  in  [4]  showed  that  waveform 
noise  related  to  optical  frequency  fluctuations  was  much  larger  for  the  <t>=7t  case  than  for  the 
0=0  case.  This  difference  was  understood  by  observing  that  largest  waveform  noise  is 
expected  at  temporal  positions  where  adjacent  input  pulses  provide  equal  contributions  and 
hence  largest  interference,  which  occurs  exactly  between  the  original  pulses  (at  T/2).  This 
time-domain  picture  provides  a  successful  qualitative  explanation  of  the  experimental 
observations. 

In  this  paper,  we  present  the  first  quantitative  investigation  of  the  impact  of  optical  comb 
frequency  offsets  on  the  time-domain  intensity  waveforms  generated  via  spectral  line-by-line 
pulse  shaping.  Instead  of  a  harmonically  mode-locked  laser,  we  employ  an  optical  comb 
generated  by  a  phase-modulated  CW  (PMCW)  laser.  By  shifting  the  center  wavelength  of  the 
CW  laser,  the  wavelength  of  the  entire  comb  can  be  shifted  accordingly,  while  leaving  the 
relative  powers,  spacing  and  phases  of  the  lines  unchanged.  PMCW  combs  are  simple  in 
configuration,  stable  even  at  the  high  repetition  rates  of  interest  here,  and  provide  much  more 
flexible  tuning  over  repetition  rate  and  comb  offset  frequency  compared  to  harmonically 
mode-locked  lasers.  Broad  PMCW  combs  have  seen  significant  recent  progress  [5].  Further, 
line-byline  shaping  on  a  PMCW  comb  has  been  reported  previously  [6],  demonstrating  the 
comb  stability  necessary  for  optical  waveform  processing.  Together  these  attributes  of  the 
PMCW  laser  enable  quantitative  study  of  the  linkage  between  variations  in  comb  frequency 
and  pulse  shaper  output. 

Fig.  1A  shows  the  schematic  experimental  setup.  A  tunable  CW  laser  is  phase-modulated 
using  a  lithium  niobate  phase  modulator  (12.5  GHz  bandwidth,  V„  ~5V  at  10  GHz)  driven  at 
frcp=9  GHz  (~1.6  V„  peak  to  peak)  to  obtain  an  optical  comb  with  lines  spaced  by  9  GFIz  as 
shown  in  Fig.  IB.  The  frequencies  are  defined  relative  to  the  CW  laser  frequency.  The  comb 
spectrum  is  observed  using  an  optical  spectrum  analyzer  (OSA)  with  0.01  nm  resolution. 
Lines  at  frequencies  higher  than  the  CW  laser  are  all  in  phase  as  measured  using  the  method 
reported  in  [7],  An  expanded  view  of  a  few  lines  is  shown  in  the  inset  on  a  linear  scale  and 
normalized  to  line  {1}.  Using  the  pulse  shaper  to  select  the  two  lines  marked  {1,  2}  and 
detuning  the  center  wavelength  of  the  CW  laser  while  the  modulation  frequency  frcp  is  fixed, 
we  are  able  to  emulate  the  effects  of  the  frequency  offsets  in  the  mode-locked  laser  of  [4]. 
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Fig.  1  (A)  Schematic  of  experimental  setup.  PM:  phase  modulator;  PS:  reflective  line-by-line 
pulse  shaper;  DET:  to  OSA  or  sampling  scope.  (B)  Optical  spectrum  of  the  phase-modulated 
CW  laser  frequency  scale  relative  to  the  CW  line.  The  frequency  lines  are  separated  by  9GHz. 

Inset:  lines  circled  are  expanded  and  plotted  on  a  linear  scale.  Lines  jl,  2}  are  selected  by  the 
line-by-line  shaper. 

Fig.  2  shows  experimental  results  where  only  two  lines,  {1,  2},  are  transmitted  through  the 
line-by-line  shaper.  Fig  2(A,C)  show  the  spectra  obtained  by  detuning  the  PMCW  comb  with 
frequency  offset  values  of  {0,  7,  14,  21, 28,  35, 42,  and  49}%  with  either  0=0  or  0=ji  applied 
to  line  {1},  respectively.  The  arrows  in  Fig.  2  indicate  the  initial  line  positions  (zero  offset) 
and  the  direction  of  detuning  as  the  frequency  offset  is  applied.  Fig.  2(B,D)  show  time-domain 
intensity  waveforms  corresponding  to  the  frequency  offsets.  Intensity  waveforms  are  obtained 
by  sending  shaped  spectral  lines  to  a  60  GHz  photo-detector  whose  output  is  measured  by  a 
sampling  scope  with  50  GHz  bandwidth.  The  vertical  arrows  indicate  the  evolution  of  the 
intensity  as  the  optical  frequency  offset  increases.  Throughout  this  paper,  intensity  waveforms 
are  normalized  to  the  peak  value  measured  with  zero  frequency  offset.  The  temporal  delay  of 
intensity  waveforms  comprised  of  two  frequency  lines  as  a  function  of  the  relative  phase 
setting  is  related  through: 

r(d>)  =  ^T  (1) 

2  n 

where  T=  l/frcp  is  the  waveform  periodicity  and  ®  is  the  phase  applied  to  one  line.  x(O)  due  to 
the  phase  settings  (0=0, tt)  are  shown  with  dashed  lines 

With  the  same  value  of  optical  frequency  offset,  peak  intensities  of  the  waveforms  of  Fig. 
2D  (®=7c)  change  much  more  strongly  with  optical  frequency  variations  than  those  of  Fig.  2C 
(0=0).  This  finding  is  consistent  with  the  results  obtained  in  [4]  using  a  non-stabilized 
harmonically  mode-locked  laser.  The  implication  is  that  the  phase  control  in  pulse  shaping  is 
the  primary  reason  for  this  difference.  In  addition,  dramatic  changes  in  waveform  shape  are 
observed  in  the  ®=7t  case  which  do  not  occur  for  0=0.  In  this  paper  we  will  develop  a 
frequency-domain  model  that  allows  quantitative  prediction  of  the  results  of  the  experiments 
shown  in  Fig.  2,  as  well  as  a  series  of  related  experiments. 

It  is  worth  noting  that  relatively  large  offsets  are  intentionally  introduced  in  this  paper,  so 
that  the  impact  of  frequency  instability  on  time-domain  waveform  generation  can  easily  be 
observed  with  a  sampling  scope.  However,  similar  time-domain  noise  effects  are  expected  to 
be  observable  even  for  smaller  frequency  fluctuations  (offsets),  provided  that  more  sensitive 
diagnostic  tools  are  available. 
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Fig.  2  Experimental  results  by  detuning  the  PMCW  comb  with  relative  frequency  offsets  of  {0, 
7,  14,  21,  28,  35,  42  and  49}%.  (A,  C)  Spectra  (linear  scale  in  relative  frequency)  for  0=0  and 
0=71,  respectively.  The  horizontal  arrows  indicate  the  initial  line  positions  and  the  direction  of 
detuning.  (B,  D)  Time-domain  intensities  for  corresponding  frequency  offsets  for  0=0  and 
0=71,  respectively.  The  vertical  arrows  indicate  intensity  peak  variations  as  the  optical 
frequency  offset  increases.  Zero  delay  positions  are  marked  by  dashed  lines.  (Video  file  of  0.96 
Mb  for  0=0  and  0.98  Mb  for  0=:r.) 


The  rest  of  this  paper  is  organized  as  follows.  In  section  2,  a  rigorous  frequency-domain 
model  is  presented  and  numerical  fit  to  experimental  data  is  provided.  Moreover,  correlations 
between  time  and  frequency-domain  pictures  are  justified  experimentally.  In  section  3,  time- 
domain  waveforms  that  show  enhanced  sensitivity  for  monitoring  of  frequency  offset 
fluctuations  are  designed  and  experimentally  demonstrated. 


2.  Modeling  and  fitting 

In  the  examples  shown  in  Fig.  2,  it  is  apparent  that  the  observed  time-domain  waveform 
variations  are  linked  to  phase  manipulation  in  the  pulse  shaping  process.  We  now  present  a 
model  which  explains  these  observations  in  the  frequency-domain. 


a.  Modeling 

For  a  grating  based  free-space  pulse  shaper,  the  effective  frequency-domain  filter  function 
H(co)  that  characterizes  its  response  as  a  linear  filter  can  be  obtained  by  convolving  the 
programmable  spatial  masking  function  M(x)  defined  by  the  liquid  crystal  modulator  (LCM) 
with  the  Gaussian  intensity  profile  of  a  single  optical  frequency  beam: 

-2(.v-g<ii); 

H{co)  =  (~/2  f M(x)e  ““  dx  (2) 

nw0  J 

where  w0  is  the  Gaussian  beam  radius  (half-width  at  1/e2  of  intensity)  and  a  is  the  spatial 
dispersion  with  units  cm  (rad/s)'1  [2].  In  the  experiment  the  frequency  lines  have  a  9  GHz 
spacing,  with  each  line  having  a  Gaussian  spatial  profile  of  radius  w(l=75  pm.  Each  LCM 
pixel  has  a  100  pm  width  and  a  corresponding  frequency  span  of  4.5  GHz;  therefore 
frequency  lines  are  centered  on  every  second  LCM  pixel.  Our  line-by-line  shaper  has 
a=3.54e-13  cm  (rad/s)'1  and  a  full-width  at  half  maximum  spectral  resolution  of  5.31GHz.  The 
spectral  resolution  is  limited  here  by  the  width  of  turning  on  one  single  LCM  pixel.  This 
obtained  resolution  ensures  observation  of  temporal  overlapping  contributions  from  adjacent 


pulses.  The  output  fields  of  the  linear  filter  Eoul(w)  is  the  product  of  the  input  field  Ej „(w)  and 
the  filter  function  H(o): 


Eou,  (a>)  =  Ejn(co)H  {(o)  (3) 

Finally,  the  time-domain  envelope  function  eoul(t)  is  obtained  by  Fourier  transforming  the 
output  field: 

\Eou,  (, o)ejmdco  (4) 

2n  J 

In  Fig.  3A,  a  schematic  representation  of  the  optical  comb  lines  in  absolute  position  (black 
solid)  and  with  an  initial  frequency  offset  £0  (blue  dash)  are  both  depicted  as  arrows.  The 
comb  line  positions  can  be  described  as  [1]: 

fm  =  mfrcp  +  S  (5) 

where  frcp  is  the  repetition  frequency,  m  is  a  large  integer,  and  e  is  the  absolute  frequency 
comb  offset.  For  our  purpose,  however,  we  note  that  the  pulse  shaper  acts  as  a  frequency  filter 
(schematically  shown  in  red  solid),  whose  center  frequency  provides  a  frequency  reference, 
which  we  denote  as  fs(red  dash).  To  optimally  select  two  comb  lines  via  pulse  shaping,  the 
comb  frequencies  {fm,  fm+i}  are  aligned  symmetrically  about  fs  by  fine  tuning  the  LCM 
position.  This  means  that  line  frequencies  referenced  to  fs  are: 

fn,  =  0«")/„,  (6) 

regardless  of  the  actual  offset  frequency  £0.  Fig.  3B  shows  the  two  lines  (blue  dash)  with  the 
filter  function  magnified  and  plotted  against  f  =  f-fs.  After  alignment,  if  the  comb  offset 
varies  from  its  initial  value  to  e’=e+  8s  (green  solid),  what  is  observed  by  the  pulse  shaper  is 
only  the  frequency  offset  change  8e.  Nonzero  values  of  8s  will  eventually  induce  power  loss 
as  the  lines  approach  the  edges  of  the  filter  function,  and  this  will  transform  into  waveform 
variations  in  the  time-domain.  It  is  important  to  emphasize^rhat  only  the  relative  frequency 
offset  is  important  in  line-by-line  pulse  shaping,  not  the  absolute  offset. 

In  Fig.  3C  we  show  M(x)  defined  as  a  four-pixel  passband  filter  (blue  solid);  4  pixels  are 
turned  on  with  constant  phase  across  the  aperture.  Recall  that  4  pixels  in  our  setup  correspond 
to  2  comb  lines.  The  effect  of  convolving  a  finite  Gaussian  beam  (green  dash)  with  a 
rectangular-like  mask  function  will  smear  out  the  abrupt  amplitude  transitions  at  the  pixel 
edges  resulting  in  a  non-rectangular  |H(co)|,  as  shown  in  Fig.  3D.  Filter  functions  on  a 
logarithmic  vertical  scale  obtained  experimentally  by  sweeping  the  wavelength  of  a  CW  laser 
(blue  solid)  and  by  simulation  (green  dash)  are  shown  for  comparison.  The  two  selected 
spectral  lines  are  indicated  by  red  arrows. 
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Fig.  3  (A)  Optical  comb  in  absolute  frequency  position  (black  solid  arrows)  and  having  an 
offset  e  (blue  dashed  arrows).  The  schematic  filter  function  is  shown  in  red.  (B)  Two  spectral 
lines  selected  by  the  pulse  shaper  are  plotted  against  relative  frequency  after  aligning  the  LCM 
(no  additional  offset:  blue  dash;  with  additional  offset  of  5e:  green  solid).  Mask  function  (solid 
line)  and  Gaussian  intensity  profile  with  radius  w„=75um  (dashed  line)  for  0=0  (C)  for  both 
lines  and  for  0=n  (E)  to  one  line.  Resulting  4-pixel  passband  filter  function  |H(m)|!:  (D) 
experimental  (solid  line)  result  by  sweeping  the  wavelength  of  a  CW  laser  and  simulated 
(dashed  line)  for  0=0;  (F)  experimental  (solid  line)  and  simulated  (dashed  line)  for  0=ji.  The 
arrows  schematically  indicate  the  comb  line  positions. 

In  Fig.  3E,  the  lower  frequency  line  (two  left  pixels)  has  an  applied  n  phase  shift,  resulting  in 
a  negative  mask  field.  3F  shows  the  resulting  experimental  (blue  solid)  and  simulated  (green 
dash)  |H(co)|  spectra  in  log  scale  -  both  reveal  an  obvious  spectral  dip.  This  dip  can  be 
understood  as  the  field  interference  due  to  an  abrupt  spatial  (7t-0)  phase  transition;  this  results 
in  diffraction  of  energy  out  of  the  main  beam,  which  is  eventually  eliminated  while  coupling 
back  into  a  single-mode  optical  fiber  at  the  output  of  the  pulse  shaper. 

In  Fig.  3F,  the  applied  phase  control  induces  a  frequency-dependent  transmission  dip  at 
the  center  of  the  passband.  This  dip  functions  as  a  frequency  discriminator  that  is  not  present 
in  the  four-pixel  passband  filter  with  constant  phase.  This  is  the  essential  point  used  to  explain 
the  observed  time-domain  waveform  instability  shown  in  Fig.  2D  when  spectral  phase  control 
is  applied.  Under  line-by-line  phase  control,  H(co)  is  no  longer  flat  in  the  frequency-domain; 
thus  a  relative  frequency  offset  er  in  the  optical  comb  will  effectively  be  transferred  to  an  extra 
energy  loss  by  the  central  frequency  discriminator.  This  loss  viewed  in  the  frequency-domain 
will  result  in  large  peak  intensity  variation  in  the  time  domain. 


b.  Data  fitting  and  discussions 


In  Fig.  4,  a  fit  is  shown  to  the  experimental  data  presented  in  Fig.  2  using  the  frequency- 
domain  model  presented  above.  Fig.  4(A,C)  show  the  normalized  optical  spectral  line 
amplitudes  in  linear  scale  with  experimental  frequency  offset  values.  The  4-pixel  bandwidth 
filter  functions  are  included  (shown  in  red  dash)  for  rt>=0  (A)  and  <t>=7i  (C)  applied  to  line  { 1 }. 
Relative  line  amplitudes  of  the  0,  1,  2  and  3  lines  are  determined  from  the  spectrum  shown  in 
Fig.  IB.  The  necessity  of  using  4  comb  lines  is  to  reflect  the  fact  that  although  only  lines  {1,2) 
are  selected  at  zero  offset,  additional  lines  from  either  side  of  the  comb  may  participate  in 
shaping  for  a  large  frequency  offset  (lines  {0  or  3},  depending  on  the  direction  of  frequency 
offset).  Arrows  indicate  the  initial  line  positions  and  the  direction  of  increasing  frequency 
offsets.  Fig.  4(B,D)  show  time-domain  intensity  waveforms  corresponding  to  the  various 
frequency  offsets  for  0=0  (B)  and  (D).  The  arrows  indicate  the  evolution  of  amplitudes 
according  to  the  change  in  frequency  offset.  The  temporal  delays  of  the  waveforms  (t)  as  a 
result  of  the  phase  manipulation  are  also  shown  (dashed  lines). 

Comparing  theoretical  calculations  shown  in  Fig.4  and  experimental  results  shown  in 
Fig. 2  demonstrate  excellent  agreement.  Several  key  features  for  discussion  are:  (i)  For  the 
0=0  case  (Fig.  2A  and  Fig.  4A),  only  the  amplitude  of  line  {1}  is  attenuated.  The  amplitude 
of  line  {2}  is  not  attenuated  due  to  a  change  in  frequency  offset  since  the  filter  function  is  flat. 
The  corresponding  time-domain  intensities  (Fig.  2B  and  Fig.  4B)  are  therefore  relatively 
insensitive  to  frequency  offsets.  Less  than  35%  variation  in  peak  intensity'  occurs  even  for  a 
49%  frequency  offset,  (ii)  For  the  0=7t  case  (Fig.  2C  and  Fig.  4C,  t(7t)=T/2),  the  amplitudes 
of  both  lines  {1  and  2}  are  attenuated  monotonically  due  to  the  increasing  frequency  offsets, 
due  to  the  central  transmission  dip.  The  corresponding  time-domain  intensities  (Fig.  2D  and 
Fig.  4D)  are  therefore  much  more  sensitive  to  change  in  frequency  offset.  For  a  49% 
frequency  offset,  the  change  in  I(t(®))  is  in  excess  of  95%.  These  two  examples  demonstrate 
the  utility  of  the  frequency-domain  model,  and  provide  a  link  with  the  time-domain  picture 
provided  in  [4].  In  particular,  the  time-domain  picture  predicts  maximum  intensity  variation 
happens  at  odd  multiples  of  T/2  (the  temporal  position  centered  between  the  original  pulses), 
where  contributions  from  adjacent  pulses  overlap  and  interact  most  strongly.  With  a  phase 
step  of  0=rc  between  adjacent  spectral  lines,  the  peak  intensities  are  delayed  to  exactly  odd 
multiples  of  T/2,  therefore  producing  a  large  intensity  variation. 
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Fig.  4  Fit  to  Fig.  2  using  the  frequency-domain  model.  (A,  C)  Spectra  fit  (linear  scale)  for  8c  of 
{0,  7,  14,  21,  28,  35,  42  and  49}%.  The  arrows  indicate  the  initial  line  positions  and  the 
direction  of  detuning.  Experimental  time-domain  traces  for  corresponding  frequency  offsets  are 
shown  for  0=0  (B)  and  0=ti  (D).  The  arrows  indicate  intensity  peak  variations  associated  with 
these  frequency  offsets.  Delay  positions  are  marked  by  dashed  lines.  (Video  file  of  0.97  Mb  for 
0=0  and  1 .00  Mb  for  O-tt.) 


Examination  of  the  0=71  amplitude  variation  (Fig.  2D  and  Fig  4D)  shows  that  the 
periodicity  of  the  waveforms  gradually  evolves  to  T/2  as  the  frequency  offset  increases.  The 


origin  of  this  repetition  rate  doubling  is  explored  in  Fig.  5.  In  Fig.  5A,  comb  lines  (arrows) 
and  the  filter  function  (red  dash)  are  schematically  sketched  against  relative  frequency.  The 
blue  solid  arrows  show  the  initial  comb  line  positions  and  the  green  dashed  arrows  show  the 
comb  with  a  frequency  offset  8s,  while  keeping  line  amplitudes  fixed.  The  lines  are  spaced 
apart  by  frcp.  When  one  of  the  spectral  lines  has  an  applied  tt  phase  shift  relative  to  the  other 
line,  the  field  of  the  filtered  comb  is  depicted  in  Fig.  5B,  with  one  line  having  negative 
amplitude.  The  corresponding  time-domain  intensity  is  shown  as  the  blue  solid  trace  in  Fig. 
5D,  a  cosine  function  having  period  of  T  with  delay  i(ct>)=T/2.  Fig.  5C  shows  the  resulting 
filtered  comb  with  offset  8e,  constituting  three  lines  denoted  as  {a,  b,  c}.  Note  that  the 
amplitudes  of  all  the  lines  are  affected  by  the  frequency-dependent  transmission  of  the  filter, 
while  the  signs  read  {+,+,-},  respectively.  The  corresponding  time-domain  intensity  for  the 
lines  with  this  particular  phase  sequence,  shown  as  the  green  dashed  trace  in  Fig.  5D,  yields  a 
shifted  cosine  function  now  with  periodicity  ofT/2.  This  phenomenon  can  be  understood  from 
the  beating  between  each  two  lines:  the  9  GHz  beating  of  lines  {a,  b}  almost  cancels  with  that 
from  {b,  c};  what  is  left  is  the  remaining  18  GHz  beating  between  lines  {a,  c}.  The  ability  to 
predict  subtle  details  of  the  experiment,  such  as  the  evolution  of  the  intensity  periodicity  from 
T  to  T/2,  again  demonstrates  the  power  of  our  frequency-domain  treatment. 


Fig.  5  Schematic  of  the  concept  of  periodicity  conversion  due  to  shifting  of  the  optical 
frequency  comb  when  a  phase  shift  is  applied  to  one  comb  line.  (A)  Initial  and  the  shifted 
frequency  comb  with  fixed  line  amplitudes  depicted  by  4  arrows  (solid:  initial;  dashed:  shifted). 
The  filter  function  is  shown  in  red  dashed  line,  (B)  Filtered  comb  without  frequency  shift.  (C) 
Filtered  comb  with  a  large  frequency  shift.  (D)  Simulated  time-domain  results  showing 
periodicity  conversion:  from  initial  comb  without  frequency  shift  (solid)  and  the  shifted  comb 
(dash). 


c.  Correlating  frequency  and  time-domain  pictures:  intermediate  phases 
& 

Thus  far  phase  steps  between  adjacent  optical  frequency  comb  lines  have  been  limited  to  <t>=0 
and  n.  This  means  the  intensity  waveform  is  shifted  from  0  (temporal  position  with  least 
impact  to  comb  offset)  to  T/2  (position  which  is  most  vulnerable  to  offset)  according  to 
temporal  delay  x(<t>)  defined  in  Eqn.  (1).  In  this  section  we  examine  other  phase  steps  in  order 
to  investigate  the  dependence  of  the  maximum  peak  amplitude  variation  on  both  frequency 
offset  and  delay. 

Experimental  and  calculation  results  together  with  <t>=0  and  n  cases  demonstrated 
previously  are  shown  in  Fig.  6.  Fig.  6A  shows  experimental  (symbols),  and  calculated  (lines), 
intensity  values  I(x(<t>))  at  5s  of  {0,  7,  14,  21,  28,  35,  42  and  49}%.  Phase  controls  are 
applied  between  0  and  n  in  rc/4  increments  to  line  {1}.  The  magnitudes  of  intensity  variation 
increase  with  increasing  t(®)  (in  the  direction  of  the  black  vertical  arrow).  Fig.  6B  shows 
l(x(0))  at  5s=49%  for  different  <t>  values  using  a  PMCW  comb  (blue  circle:  experimental; 
blue  solid  line:  calculation).  x(O)  is  also  included  within  the  figure:  experimental  (diamonds) 


and  calculation  (dotted  line).  The  perfect  linearity  of  the  trace  indicates  excellent  control  over 
phase  settings  to  individual  lines.  In  all  calculations,  the  optical  comb  line  amplitudes  follow 
the  experimental  spectrum  as  shown  in  Fig.  IB.  Results  show  errors  of  less  than  10%  between 
experimental  data  and  calculations.  Together,  Fig.  6A  and  6B  indicates  that  the  amplitude 
change  of  l(x(0))  increases  monotonically  both  with  frequency  offset  and  temporal  delay,  in 
agreement  with  the  time-domain  view  presented  previously. 

As  a  further  example,  calculated  results  of  l(t(<t>))  at  8e=49%  assuming  an  ideal  comb 
comprised  of  equal  amplitude  lines  (red  square)  are  plotted  against  cos2(<t>/2)  (red  dashed  line) 
in  Fig.  6B.  Using  the  linear  dependence  of  t(O)  on  <J>,  we  infer  that  for  an  ideal  comb,  the 
range  of  intensity  variation  with  frequency  offset  follows  a  temporal  mask  of  the  form 
cos2(jtfrcpT).  We  will  use  this  empirical  observation  as  an  aid  to  interpret  some  of  our  later 
results. 

Please  note  that  the  form  of  this  relation  may  depend  on  pulse  shaper  resolution.  In  the 
current  experiments,  our  pulse  shaper  spectral  resolution  is  1.7  times  finer  than  the  line 
spacing,  which  leads  to  a  regime  where  the  waveform  interference  in  the  time-domain  is 
dominated  by  contributions  from  only  nearest  neighbor  pulses.  This  corresponds  most  closely 
to  the  simple  time-domain  model  put  forth  in  [4], 
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Fig.  6  (A)  Experimental  (symbols:  <t>=  {0, 7t/4,  n/2,  3tt/4  and  n})  and  calculated  (lines:  d>=  0  to 
7t  in  7t/4  increments)  waveform  amplitudes  at  corresponding  t(O)  for  6a  of  (0,  7,  14,  21 , 28,  35, 
42  and  49}%.  (B)  Intensity  values  for  49%  frequency  offset  at  corresponding  t(<!>)  for  different 
<!>  values  using  phase-modulated  CW  comb  (blue  circle:  experimental;  blue  solid  line: 
calculation).  Calculated  (dotted  line)  and  experimental  (diamonds)  t(O)  is  also  included. 
Intensity  values  assuming  ideal  comb  lines  (red  square)  is  plotted  with  cos!(<I>/2)  (red  dashed 
line). 


3.  Tirrfe-domain  waveforms  designed  for  monitoring  optical  frequency  fluctuations 

a.  Design  approach 

As  discussed  above,  phase  manipulation  of  one  spectral  line  for  a  two-line  cosine  time-domain 
waveform  is  sensitive  to  spectral  frequency  shifts.  We  can  therefore  ask  the  question:  can 
line-by-line  pulse  shaping  be  used  to  monitor  variations  in  the  comb  offset  frequency,  at  least 
in  the  large  fluctuation  regime?  For  the  purpose  of  frequency  fluctuation  monitoring,  it  would 
be  useful  to  design  waveforms  that  not  only  exhibit  temporal  positions  that  are  sensitive  to 
frequency  fluctuation,  but  also  have  stable  reference  points.  The  existence  of  a  reference  point 
that  is  insensitive  to  optical  frequency  fluctuations  would  allow  one  to  distinguish  waveform 
changes  arising  from  pure  frequency  fluctuations  from  the  effects  of  laser  intensity 
fluctuations.  This  should  also  increase  the  sensitivity  of  such  a  scheme  for  monitoring  optical 
comb  frequency  fluctuations. 

Based  on  the  previous  discussions,  one  simple  choice  of  target  waveform  would  be  a  time- 
domain  cosine  generated  by  selecting  two  spectral  comb  lines.  In  general,  one  can  select  any 


two  spectral  lines  separated  by  Nfrcp,  where  N  is  an  integer.  For  the  N=1  example  discussed 
above,  we  have  shown  that  there  is  significant  sensitivity  to  frequency  fluctuations  at  time 
positions  equal  to  odd  multiples  of  T/2  but  not  at  multiples  of  T.  Let  us  assume  that  this 
behavior  applies  also  for  N>1,  which  is  justified  by  the  results  given  below.  Then  any  two 
lines  spaced  apart  by  2frL.p  or  greater  will  exhibit  waveform  positions  both  sensitive  and 
insensitive  to  comb  fluctuations,  since  there  will  be  peaks  situated  at  multiples  of  T/N. 
Furthermore,  applying  a  phase  shift  onto  one  of  the  two  lines  allows  us  to  finely  shift  the 
positions  of  the  peaks  for  optimum  sensitivity. 

To  illustrate  this  argument,  Fig.  7  shows  simulation  results  using  the  frequency-domain 
model.  For  simplicity,  an  ideal  comb  of  fr<;p=9  GHz  spacing  (T=l  11.1  ps)  is  assumed  with 
spectral  comb  lines  of  equal  amplitudes.  Fig.  7A  shows  a  filter  designed  to  select  two  lines 
spaced  by  2frep  (N=2)  which  removes  the  intermediate  line.  Optical  frequency  comb  offsets  of 
{0-50}%  with  10%  per  step  are  indicated.  The  corresponding  time-domain  waveforms  are 
shown  in  Fig.  7B.  1(0)  is  invariant  while  i(T/2)  is  sensitive  to  frequency  fluctuations.  The 
cos2(7ifrcpt)  is  plotted  (green  dashed  line)  to  emphasize  the  temporal  intensity  variation  relation 
discussed  in  the  previous  section.  The  ratio  of  I(T/2)  to  1(0)  offers  a  clear  indication  of 
frequency  offsets,  especially  for  large  offsets.  Fig.  7C  shows  the  semi-log  plot  of  the  contrast 
ratio  with  various  offset  values.  The  scenario  is  different  if  one  of  the  lines  is  71-shifted,  shown 
in  Fig.  7D.  Since  the  waveforms  are  delayed  by  T/4,  all  the  peaks  encounter  the  same  amount 
of  intensity  variation.  One  interesting  observation  is  how  the  cos2(7tfrcpt)  function  nicely  cuts 
the  intensity  waveform  at  50%  offset.  This  once  again  supports  our  conclusion  drawn  from  the 
previous  section. 
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Fig.  7  Two  lines  selected  at  2frcp  by  blocking  the  center  line.  (A)  Ideal  comb  lines  shifted  with 
5e  { 0-50} %  in  10%  increment  with  no  phase  control.  (B)  Resulting  time-domain  intensity: 
peak  at  time  zero  is  stable  while  peak  at  time  T/2  is  sensitive  to  offsets.  (C)  Semi-log  plot  of 
I(T/2)/I(0).  (D)  Ideal  comb  lines  shifted  with  offsets  {0-50}%  in  10%  increment  with  <t>=7r  to 
one  line.  In  (B,  D),  green  dashed  lines  depict  cos2(7tfrcpt). 


b.  Experimental  results 

A  broader  and  flatter  optical  comb  is  favored  for  this  experiment  to  better  approach  the  ideal 
comb  (equal  intensities  for  the  comb  lines).  A  low  V„  lithium  niobate  phase  modulator 
(20GHz  bandwidth,  V„  of  -2.8V  at  1  GHz)  is  modulated  at  frcp=9GHz  with  a  driving  voltage 
of -2.6  V„  peak  to  peak  to  obtain  the  spectrum  shown  in  Fig.  8A.  Relative  phases  of  the  lines 
are  measured  for  lines  {-6  to  5}.  Lines  {1,  0,  -3,  and  -5}  are  k  out  of  phase  with  the  rest, 
which  agree  with  calculation  using  a  Bessel  series  [8].  Fig.  8B  shows  an  expanded  portion  of 
the  comb  on  a  linear  scale  with  line  powers  normalized  to  line  {3}.  Lines  {2  and  4}  are 
selected  to  pass  through  the  pulse  shaper  while  the  others  are  blocked  by  applying  amplitude 
control  to  the  line-by-line  shaper  to  obtain  two  spectral  lines  spaced  by  2frcp. 


Frequency  (GHz)  Frequency  (GHz) 

Fig.  8  (A)  Spectrum  obtained  from  a  low  V„  phase  modulator.  Circled  spectral  portion  is 
expanded  in  (B),  in  linear  scale  and  normalized  to  line  {3}.  Lines  {2,  4}  are  selected  by  the 
line-by-line  shaper  while  the  others  are  suppressed. 

Fig.  9  shows  experimental  and  calculated  results  for  this  waveform.  Fig.  9(A,B)  show  the 
experimental  spectra  and  line  amplitudes  used  in  the  calculation  with  8e  of  {0,  7,  14,  21,  28, 
35,  42,  and  49}%,  both  in  linear  scale.  In  both  figures,  the  horizontal  arrows  indicate  the 
initial  line  positions  and  the  direction  of  detuning.  In  Fig.  9B,  the  amplitudes  of  the  spectral 
lines  reflect  the  experimental  spectrum  shown  in  Fig.  8B.  Fig.  9(C-E)  show  resulting  time- 
domain  profiles  for  the  various  shifts  of  the  frequency  comb  with  phase  ®={0,  rr/2,  and  7t} 
applied  to  line  {2},  respectively.  Fig.  9(F-H)  show  the  corresponding  calculated  intensities 
reflecting  the  experimental  conditions.  In  Fig.  9(C-H),  vertical  arrows  indicate  intensity 
variations  as  the  optical  frequency  offset  increases  for  different  temporal  positions;  zero  delay 
and  T/2  positions  are  marked  as  an  aid  to  the  eye. 

.  Again  the  agreement  between  data  and  theory  is  excellent.  Even  the  detailed  features 
observed  at  large  frequency  offsets  show  excellent  correspondence.  For  example,  both  C  and 
F  show  the  repetition  period  changing  from  T  to  T/2;  both  D  and  G  show  temporal  shifts  in 
the  peak  near  T/2;  and  both  E  and  H  show  symmetric  waveform  shifts  about  t=0.  In  the  end,  it 
is  the  waveform  with  0=0  (C  and  F)  that  is  of  primary  interest  for  monitoring  of  optical  comb 
frequency  shifts.  In  contrast  to  Fig.  7B,  where  intensity  values  at  1(0)  are  invariant  to 
frequency  offsets,  1(0)  in  Fig.  9(C,F)  increases  slightly  with  increasing  frequency  offset.  This 
results  in  an  enhanced  sensitivity  to  frequency  shifts  and  is  attributed  to  the  unevenness  of  the 
PMCW  line  powers  in  Fig.  8B.  Under  large  frequency  offset,  the  strong  central  line  {3}  that  is 
initially  blocked  contributes  to  pulse  shaping  and  accelerates  the  variation  in  the  intensity 
waveforms. 
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Fig.  9  (A,  B)  Experimental  spectra  and  the  lines  used  for  calculation  (linear  scale)  for  two  lines 
with  spacing  of  2frep  with  8c  of  (0,  7,  14,  21,  28,  35,  42,  and  49)%.  The  horizontal  arrows 
indicate  the  initial  line  positions  and  the  direction  of  detuning.  (C-E)  Experimental  time- 
domain  waveforms  for  <t>=  {0,  rt/2,  and  jt),  respectively.  (F-H)  Calculated  time-domain 
waveforms  for  <t>=(0,  7t/2,  and  it),  respectively.  The  vertical  arrows  indicate  intensity 
variations  as  the  optical  frequency  offset  increases  for  different  temporal  positions.  Temporal 
positions:  dash-dot:  zero  delay;  dashed:  T/2. 


c.  Larger  N  values 

One  may  extend  the  spacing  between  the  two  selected  lines  to  achieve  similar  results.  Fig.  10 
shows  simulation  results  for  an  ideal  (equal  amplitude)  comb  for  N=3  (A-D)  and  N=4  (E-FI). 
Note  how  well  the  frequency-offset  intensities  correspond  to  the  cos2(7tfrcpt)  function.  The 
time-domain  picture  is  clearly  illustrated:  frequency  sensitivity  of  waveforms  generated  via 
line-by-line  pulse  shaping  in  the  100%  duty  cycle  regime  is  least  noticeable  at  time  positions 
equal  Jto  0  (or  integral  multiple  of  the  repetition  period  T)  and  increasingly  noticeable  for 
increasing  time  offsets.  Maximum  intensity  variations  with  frequency  offsets  occur  for  time 
positions  of  T/2  and  its  odd  multiples. 
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Fig.  10  Simulation  spectra  for  two  lines  with  N=3  (A)  and  N=4  (E)  with  5e  of  (0-50)%  in  10% 
increments.  Resulting  time-domain  waveforms  are  given  for  N=3  (B-D)  and  N=4  (F-H)  with 
phase  shifts  applied  to  one  line  <t>={0,  n/2,  n},  respectively.  Green  dashed  lines  depict 
cos'(jtfrt.nt).  Temporal  positions:  dash-dot:  zero  delay;  dashed:  T/2. 


5.  Conclusion 

In  summary,  the  impact  of  frequency  comb  shifts  on  time-domain  waveforms  generated 
through  line-by-line  shaping  has  been  quantitatively  investigated.  The  influence  of  frequency 
offset  on  generated  waveforms  is  explained  from  both  time  and  frequency-domain 
perspectives.  A  quantitative  frequency-domain  model  is  provided  starting  from  pulse  shaper 
theory.  Time-domain  intensity  variations  that  have  previously  been  observed  in  pulse  shaping 
due  to  an  unstable  mode-locked  frequency  comb  are  experimentally  emulated  using  a 
wavelength  detuned  phase-modulated  CW  laser,  which  allows  quantitative  analysis.  The 
experimental  data  are  in  excellent  agreement  with  the  predictions  of  our  frequency-domain 
model.  Finally,  we  propose  and  model  pulse  shaper  waveforms  that  should  be  useful  for 
monitdring  the  optical  frequency  fluctuations  of  the  input  comb. 
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